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Abstract 

Motivated by ongoing developments at HI7S at TUNL that include increased photon flux and the 
ability to circularly polarize photons, we calculate several beam-polarization/target-spin dependent 
observables for elastic Compton scattering on the deuteron. This is done at energies of the order of 
the pion mass within the framework of Heavy Baryon Chiral Perturbation Theory. Our calculation 
is complete to 0{Q^) and at this order there are no free parameters. Consequently, the results 
reported here are predictions of the theory. We discuss paths that may lead to the extraction 
of neutron polarizabilities. We find that the photon/beam polarization asymmetry is not a good 
observable for the purpose of extracting q„ and (3n- However, one of the double polarization 
asymmetries, T,^, shows appreciable sensitivity to 7i„ and could be instrumental in pinning down 
the neutron spin polarizabilities. 
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I. INTRODUCTION 



Compton scattering studies on the nucleon and deuteron have shed hght on the low- 
energy dynamics of quarks inside a nucleon. An electromagnetic probe by nature, Compton 
scattering captures information about the response of the charge and current distributions 
inside a nucleon to a quasi-static electromagnetic field. These responses are quantified in 
terms of the nucleon polarizabilities. 

To lowest order in photon energy, the spin-averaged amplitude for Compton scattering 
on the nucleon is given by the Thomson term 

Amp = — —t ■ e', (1) 

where Ze, M represent the nucleon charge and mass respectively and e, e' specify the po- 
larization vectors of the initial and final photons respectively. At the next order in photon 
energy, contributions arise from electric and magnetic polarizabilities — ole and 13m — which 
measure the response of the nucleon to the application of quasi-static electric and magnetic 
fields. The spin-averaged amplitude is expressed as: 

Amp = e ■ e' f — — + ujuj' Aixole] + e x k ■ e' x k' ujuj'AtiPm + 0{uj'^) . (2) 

Here, k^ = {uj,k), k'^ = {uj',k') specify the four-momenta of the initial and final photons 
respectively. The associated differential scattering cross section on the proton is given by 
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da _ / y fuj^ 



-(1 + cos^i 



— - — i^-{aE + Pm){i + COS ey + -{aE - Pm){i - COS ey ] + 



(3) 



It is essential to mention here that for the purpose of this paper, and /3m are the so- 
called Compton polarizabilities. In the non-relativistic limit, these reduce to to the static 
electromagnetic polarizabilities. From here on we shall drop the subscripts E and M as well 
as the bars over a and f3. 

At the next order in u, we can access the spin polarizabilities (7's) of the nucleon. Unlike 
a and /3, however, there is no simple classical definition for the spin polarizabilities. Ref. 0] 
provides a treatise on polarizabilities for introductory purposes. 

Effective Field Theory (EFT) techniques have been successfully used to explore the 
infrared regime of QCD. They are rooted in the symmetries of the underlying theory. Re- 
search programs in the recent past have facilitated the development of a framework based 
on EFT techniques that allows a consistent description of low-energy hadronic and nuclear 
processes. This framework involves building an EFT from the most general Lagrangian 
involving baryons and mesons below the chiral symmetry breaking scale, A^. In its simplest 
form, this EFT involves the lightest mesons and baryons, namely pions and nucleons, as the 
only degrees of freedom. The heavier mesons and baryons are integrated out of the theory. 
A perturbative expansion is then built as a function of the ratio of two scales - one light 
and the other heavy - to be identified with the degrees of freedom in the theory. The price 
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one has to pay with such a choice is that the domain of vahdity of the theory is hmited by 
the hghtest mass scale not exphcitly included in the theory. In this case, it is the energy 
required to excite the A(1232) isobar, which is of the order of 300 MeV. 

This EFT is known as Chiral Perturbation Theory (xPT). In this theory pions interact 
via derivative interactions and/or insertions of the quark mass matrix. The hierarchy of 
scales in this theory is set by the 'small' momentum, p, and the pion mass, which explicitly 
appears in the interactions, i.e. Q = {p,mT^), and A^. Thus the S'-matrix elements can be 
expressed as expansions in powers of the ratio of the two scales, 9iziEiIIhLl_ The dynamics 
at distances shorter than the wavelength of the probe or those involving mesons other than 
the pions are accounted for by local operators. In the purest form of xPT the strengths 
(in other words, the coefficients) of these operators, called "low-energy constants" (LECs) 
should be extracted from experiments. The resulting theory can then be used to predict the 
outcome of other experiments. Baryons are usually included in xPT by treating them as 
heavy objects because of the fact that the momenta involved are very small compared to 
the baryon rest masses. The resulting theory is called Heavy Baryon Chiral Perturbation 
Theory (HBxPT). A detailed discussion can be found in Refs. 0,1^. 

HBxPT has been very successful in the single-nucleon sector and its successes include 
Compton scattering on a proton. Nucleon Compton scattering has been studied in HBxPT 
in Ref. 0], where the following results for the polarizabilities were obtained to 0{Q^): 

2 2 

7iP = 7in = 0/3^2 2 = 4-4 X 10-Hm\ 
2 2 

72p = 72n = ,0/3^/2 2 = 2-2 X 10-^fm^ 

73P = 73. = 'Y^, , = 1.1 X 10-^fm^ 

= 74n = - ^^fsl , = -1-1 X 10-^ fm^ (4) 

Here, we have used qa = 1-26 for the axial coupling of the nucleon, and = 93 MeV as the 
pion decay constant. At this order, the polarizabilities are manifestations of isoscalar pion 
loops (see Fig. |2I), and so are equal for the proton and the neutron. Note that the polariz- 
abilities are predictions of HBxPT at this order. The 0{Q^) values of the polarizabilities 
are modified by higher-order corrections which include what can be quite sizeable effects 
due to the A isobar [1, In our calculations, the 0{Q^) values of the polarizabilities will 
serve as a reference point as we seek to understand how such shifts away from these values 
(Eq. (0])) will impact observables. 

There has been an avalanche of experimental data on unpolarized Cornpton scattering 
on the proton in the recent past at photon energies below 200 MeV 0, 13] ■ The proton 
electromagnetic polarizabilities are now quite well established. The current Particle Data 
Group (PDG) values for the proton are: 

ap = (12.0 ±0.6) X 10"^ fm^ 
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I3p = (1.9 ±0.5) X 10"^ fm^ 



(5) 



Theoretical efforts have matched experiments with calculations being done up to the fourth 
order in Q E3l- The HBxPT calculations give a very impressive overall description of 

the data for u, ^\ < 200 MeV. 

In contrast, due to the lack of free neutron targets, information about the neutron po- 
larizabilities has to be gathered through alternative methods. One such method involves 
scattering neutrons on lead to access the Coulomb field of the target and examining the 
cross-section as a function of energy. Currently there is much controversy over what this 
technique gives for Two experiments using this same technique obtained different results 
for an (Refs. [11 Q)- 

On = (12.6 ± 1.5 ±2.0) X lO^^fm^; (6) 
On = (0.6 ±5.0) X 10-^fml (7) 



Enik et al [TJ revisited Ref. ^ and recommended a value of between 7 x 10 fm^ and 
19 X 10"^ fm^. 



Also, quasi-free Compton scattering from the deuteron was measured at Mainz [IJ] for 
incident photon energies of = 236 — 260 MeV and one-sigma constraints on the polariz- 
abilities were reported to be: 

an = (7.6-14.0) X 10-^fm^ 

/?„ = (1.2-7.6) X lO^^fml (8) 

Note that it is necessary to pin down only one of a„ and In practice, one uses the 
dispersion sum rule that relates the sum of the polarizabilities to the total photoabsorption 
cross-section in order to extract the other. For the proton 

a+p=— P \ ^ dcu = (13.8 ± 0.4) X 10"^ fm^ (9) 

271-' Juj-fh OJ 



where, a* (a;) is the^ total photoabsorption cross-section and ujth is the pion-production 



threshold. In Ref. [15[ the sum rule for the neutron was extracted from the deuteron pho- 
toabsorption cross-section by subtracting the proton contribution obtaining, 

a„ + /5„ = (14.4 ± 0.66) x 10"^ fml (10) 

Given the discrepancy between the measurements of in Eqs. ^ and ((Zj), it is impor- 
tant to measure it by other means. Compton scattering on a nuclear target is an obvious 
candidate. However, processes with A > 1 pose a different kind of challenge because they 
require an understanding of the inter-nucleon interaction. For example, deuteron structure 
is governed by the NN interaction and hence, when analysing '-jd data, one has to include 
the effects of two-body currents over and above the single-nucleon amplitude. There 
are different schools of thought regarding this. The first school was that introduced by 
Weinberg ^ and later refined. Another formulation where pions are treated perturbatively 
was developed by Kaplan, Savage and Wise (KSW) However, the KSW scheme is 

known to break down at much lower energies than the Weinberg formulation. We use the 
Weinberg formulation for our calculations. Elastic Compton scattering on the deuteron has 
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been studied in this theory up to 0{Q^) |lO|, llll as a means to extract the spin-independent 
nucleon isoscalar polarizabihties. 

Recent experimental advances have made it possible to venture into the area of polar- 
ization observables. Technology to polarize targets and the current developments at the 
HI7S facility at TUNL jl^ have especially motivated research in this area, in particular 
the research reported in this paper. The HI7S upgrade program includes increased photon 
flux and the ability to circularly polarize photons. We believe that these efforts will make it 
possible to design experiments to access the spin polarizabihties of the nucleon. Hildebrandt 
et al have studied polarization observables for Compton scattering on the nucleon (both 
proton and neutron) with a focus on the spin polarizabihties. Their calculation involves free 
nucleons and is beneficial for understanding processes that involve quasi-free kinematics like 
•yd — > '-ynp as opposed to •yd — > yd. In this paper, we calculate elastic Compton scattering 
on the deuteron to 0{Q^) but focus on specific observables so as to construct road-maps 
to extract the neutron polarizabihties, a„, /5„ and 7i„ . . . 74^. These calculations are ex- 
ploratory, and are expected to provide a reference for polarization observables for elastic 
Compton scattering on the deuteron. As such they will serve to focus further explorations 
of these observables at 0{Q^) and beyond. 

Until now, there have been only limited efforts to measure nucleon spin polarizabihties. 
The only measurements attempted have been for the forward and backward spin polariz- 
abihties. The experimental values for the backward and the forward spin polarizabihties 
quoted later in this Section include the contribution from the t-channel 7r° pole. However, 
deuteron being an isoscalar, this graph does not contribute in the theory discussed in this 
paper. 

The proton backward spin polarizability, 7,7 = 7i + 72 + 274 has been extracted by the 
LEGS collaboration from unpolarized Compton scattering on the proton 0]. But, the 
extracted value 

7^p = (-27.1 ± 2.2) X 10-^fm^ (11) 

contradicts predictions of standard dispersion theory j2l[ US, l23| and xPT Hi The 
theoretical prediction for ^T^p from xPT is —36.7 x 10~^ fm^ |2^. The LEGS result 
is also in disagreement with the TA PS result 7-71^ — 

(-35.9 ± 2.3) X 10"^ fm^ 0. The 
most recent Mainz experiment ji^ using the Mainz 48 cm 0x 64 cm Nal detector and 
the Gottingen recoil detector SENECA in coincidence, extracted a backward polarizability 
value ranging from (-36.5 to -39.1) x 10"'^ fm'^, depending on the parameterization of the 
photomeson amplitudes. These values are consistent with the Mainz measurements using 
the LARA detector (27I . |28| , which extracted a backward polarizability value ranging from 
(-37.1 to -40.9) xlO~^ fm^, again depending on the parameterization of the photomeson 
amplitudes. These results are also in agreement with the earlier result from Mainz 
The neutron backward spin polarizability was determined to be 

y^n = (58.6 ± 4.0) x 10"^ fm^ (12) 



from quasi-free Compton scattering on the deuteron j29|]. This experiment also used the 
Mainz 48 cm 0x 64 cm Nal detector and the Gottingen recoil detector SENECA in coinci- 
dence. 

The forward spin polarizability 70 is related to the energy- weighted integrals of the differ- 
ence in the hehcity-dependent photoreaction cross-sections {0-1/2—0-3/2) through the following 
sum rule: 

70 = 71 - (72 + 274) = A ""'^'"s""'^' da;' (13) 
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where LVth is the pion-production threshold. The following results on 70 were estimated using 
the VPI-FA93 multipole analysis [sj: 

7op ~ -1.34 X 10-^fm^ (14) 
7o„ ~ -0.38 X 10"^ fm^ (15) 

Our goal in this paper is to calculate beam-polarization/deuteron-spin dependent observ- 
ables for the elastic Compton scattering process on the deuteron at energies below the pion 
mass with an emphasis on the role of the neutron polarizabilities. In a recent paper, Chen, Ji 
and Li [sj , performed calculations in pionless effective theory with a similar goal. However, 
that theory breaks down at ~ 50 MeV in contrast to the EFT used here which is applicable 
from around 50 MeV to the pion-production threshold. Knowing that the proton Thomson 
term dominates the jd scattering at low energies (tens of MeVs), it is essential to study the 
process at energies where the sensitivity to the neutron polarizabilities will be appreciable. 

In Section |n] we describe the calculation with focus on calculating the observables de- 
scribed in Section irm We discuss the strategy used for our calculations in Section HVl Then 
in Section |3 we report our results pertaining to these observables and also assess some 
theoretical ambiguities. We conclude in Section IVIl 



II. OUR CALCULATION 

A comprehensive explanation of the HBxPT calculation procedure for Compton scatter- 
ing on the deuteron is available in Ref . • The focus in this section is primarily on the steps 
that lead to the calculation of photon-polarization/deuteron-spin dependent observables. 




FIG. 1: The anatomy of the calculation. The irreducible am- 
plitude is computed in HBxPT and sewn to external deuteron 
wavefunctions to give the matrix element for Compton scattering 
on the deuteron. 

We proceed by calculating the irreducible amplitude for Compton scattering on the 
deuteron (the central blob in Fig. and then sew it to the incoming and outgoing deuteron 
wavefunctions to obtain the scattering amplitude, i.e. 

^= (^'jT^TVivl^d). (16) 

The irreducible amplitude, T^nn, consists of one-body interactions (see Figs. El and E)) and 
two-body interactions (see Fig. HJ, 

T^NN = T^NN + PjNN- (l*^) 
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The first term is the one-body amplitude in the scattering process where only one of the 
nucleons interacts with the photon and the second term is the two-body amplitude where 
both the nucleons take part in the Compton scattering. 

The one-body amplitude in the ''jN center-of-mass (CM) frame can be written as: 

T^^ = l^ie' ■€ + A2e' -ke-k' + lA^a ■ (e' x e) + lA^a ■ [k' x k) e' ■ e 
+iA5a -[{e' xk)e-k' ~{exk')e' ■k]+ iA^a ■ [{e' x k') e ■ k' - {e x k) e' ■ k]} . (18) 

where e(e') and k{k') are the unit vectors for polarization and momentum of the incoming 
(outgoing) photons respectively, a = 2S, where S is the quantum-mechanical spin operator 
for the nucleon. The six invariant amplitudes, Ai . . . Aq, which are functions of the photon 
energy, u, and Mandelstam t, are obtained by calculating the single-nucleon amphtudes to 
0{Q^) in HBxPT. The relevant Feynman graphs are shown in Figs. |21 and 121 







FIG. 2: Tree level diagrams for the one-body amplitude. 

Defining T = uj/niT^ and t = — 2T^(1 — cos0), where 6 is the center-of-mass angle between 
the incoming and outgoing photon momenta, one finds 0, 0, H^l : 



Ai 

A2 
As 



A, 
A5 
A, 



M + 87r/2 



Vl - T2 + 



2-t 



^ arctan - /i(T, t) 



M2 87r/2m, {-t)^^ 



UJ 



2M2 



[Z{Z + 2k) - {Z + Kfcos9] 



{2Z - l)gAm^ Tt 



l-t 



T 



arcsin^ T - T + 2T^ sin^ OhiT, t) 



[Z + k)'^uj 



[Z + Kfu {2Z - 1)^A^2 T 



2M2 



[/5(T,t)-2T2cos0/3(T,t)], 



Z{Z + k)uj ^{2Z-l)gAUj'^ T 



2M2 



87r2/>, (1 - 1) 87r2/>, 



■[J5(T,t)-2T2j3(T,t)], (19) 
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"-N. 





FIG. 3: Contribution to the one-body amplitude at These 
diagrams contain one pion loop. 



where 



with 



Jo 



dz arctan 



2x71 - T2 -2^ 



/3(T,t) = 



1 2(l-z)^-t(l-T^z^) 
4(1- W) _t(i_^)2' 

x(l — x)z{l — z 



h{T,t) - 

/5(T,i) = dx t dz 
Jo Jo 



dx / dz 
Jo 

da; / dz arcsin— — , 

JO o it 

il-zf . Tz 
arcsm— , 



'3 r Tz TzS 
arcsin h 



R 



(20) 



S = ^Jl - T^z"^ - t{l - zyx{l - x), R = ^Jl-t{l- zfxil-x). 



(21) 

At low energies (a; < m7r), these six invariant amplitudes can be Taylor expanded in cj as 
follows: 



A2 

As 

A, 



u 
2AP 



[Z{Z + 2k) -{Z + cose] + Af + (71 - (72 + 274) cos{e))uj^ + 0{uj^), 



{Z + Kfuj 

2M2 



+ 72C^' + 0(cu'), 



8 



Here, 



2M2 
2M2 



(22) 



, v4f and Af are contributions from the 7r°-pole graph, but are not relevant in 



the current calculations. The zeroth-order term in uj is the leading-order Thomson term (as 
in Eq. (P)) and the first-order term contains contributions from the anomalous magnetic 
moment. The u;^ term is dependent on a and /5 and the spin polarizabilities (7's), appear in 
the a;^ term. It is fairly evident that, to leading order in Q, the nucleon polarizabilities are 
manifestations of the pion cloud and hence any theory attempting to describe this property 
of nucleons should have pions as explicit degrees of freedom. 

Since the amplitude in Eq. (fTHj) is in the CM, it has to be boosted to the CM 
frame to obtain T^^^. Consequently, 



2„2 



boost 



2M2 



-[e-k'e' ■k + 2{e-pe' ■ k + e-k'e' -p)] 



(23) 



must be added to the one-body amplitude in the 'jN CM frame (Eq. ()18|) ) when we are 
computing Compton scattering from the two-nucleon system at 0{Q^) in the 'yd CM frame. 
This additional term is necessary to ensure that the vectors e and e' stay orthogonal to k 
and k ' respectively in the new frame [sT 



FIG. 4: Contributing two-body diagrams at 0{Q^). Permutations 
are not shown. 



The two-body diagrams that contribute to the Compton scattering process at 0{Q^) are 
shown in Fig. |3 The two-body amplitude can be expressed (in the jd CM frame) as: 



2N 
-yNN 



where 



da) 



[T ■ T - n 



2[cu2 -ml- {p-p' + i(/c + k')y 



+ (1^2) 



(24) 



(25) 
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.(,) _ 6 ■ 6 ■ {p- p' - lik -k'))a- ■ {p- p' + l{k-k')) 



tic) 



, - - , - - +(1^2), (26) 

[{p-p' -l{k-k')y + ml][{p-p' + l{k-k')y + ml] 

e'-{p-p' + \k) a' -eg'- {p- p' + \{k-k')) 

[^2 _ ^2 _ p*/ + 1 + k')Y][{p-p' + \{k- k')Y + ml] 

+ (1 ^ 2), (27) 

^(d) ^ g • (P - + |fe ^'-{P- P' -\ik- k'))a' ■ e' 

[^2 _ ^2 _ ^/ + ^/))2][(p'_ p-/ _ _ ^/))2 + ^2] 

+ (1 ^ 2), (28) 

^(,) _ 2e- (p-p^ + e' -{p-p' + \k) a' ■ {p - p' - \{k ~ k')) 

~ [uj^ -ml- {p- p' + lik + k')y][{p-p' - \{k - k')Y + ml] 

■ ip-p' + \{k-k')) , ^ 
X : ^ ^1 ^ + 1 ^ 2 , 29 

Here again, e,e',k,k' have their usual meaning. p{p') is the initial (final) momentum of 
the nucleon inside the deuteron, and a^(a^) is twice the spin operator of the first (second) 
nucleon. The numbering of the nucleons is arbitrary and hence (1 2) denotes the term 
when the nucleons are interchanged. 



III. THE OBSERVABLES 

In this and the subsequent sections we shall use the following convention for the co- 
ordinate system. The beam direction is defined as the z axis, the x — z plane is the 
scattering plane with the y axis being the normal to it. For linearly polarized photons, the 
polarization vectors in the initial state can be along the x or the y axis, i.e., e = x or y. 
For circularly polarized photons in the initial state, the beam helicity. A, can be ±1. The 
photons are right circularly polarized (RCP) if the beam helicity, A = +1 (e^ = — ^^) or 

left circularly polarized (LCP) if the beam helicity, A = — 1 (e_ = ^7^). 



A. Photon/Beam Polarization Asymmetry 

The photon/beam polarization asymmetry uses linearly polarized photons and an unpo- 
larized target. The photon polarization asymmetry is defined as: 

( da_\ / do- \ 

\dn)^~^ ydn)y 

where the subscript x{y) refers to the differential cross-section when the initial state of 
polarization of the photon is in (perpendicular) to the plane of scattering. It should be 
noted that the denominator of Eq. pO|l is simply the unpolarized total differential cross- 
section. In calculating (^)j, where i = {x,y) we sum the relevant matrix elements as: 

J2 \{M',e'\T,NN\i = lM)\\ (31) 

M' ,M,e'=x,y 
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Here, M(M') defines the initial (final) spin state of the deuteron and e' is the final polar- 
ization state of the photon. 



B. Double-Polarization Asymmetry 

The double-polarization asymmetry involves circularly polarized photons and a spin- 
polarized target. The expectation is that these observables can provide insight into the spin 
polarizabilities of the target. 

When the target is polarized along the beam direction (parallel or anti-parallel), the 
corresponding observable is called the parallel target polarization asymmetry and is defined 
to be: 

( da_ \ ( — \ 

Idcht + \lnhl 

Parallel (anti-parallel) arrows in the subscript symbolize target polarization parallel (anti- 
parallel) to the beam helicity. For A = ±1 the matrix elements for (^)tT ^I's summed 
as: 

\{M',e'\T,MN\e = ex.M=+l)\\ (33) 



taking the 2;-axis as the spin quantization axis. For (^)ti) we simply replace |M = +1) by 
|M = — 1) in the initial state. 

The target may be polarized along ±a; too. In this case, the observable is called the 
perpendicular polarization asymmetry and is given by: 

^x,(A=±l) - Trf^T —7d^^ • K^V 

ydo,'^-^ ^dnn*- 

Again, the direction of the second arrow in the subscript denotes the target polarization 
is along the +x or —x directions. We have to be careful in defining the matrix elements 
because the spin state of the deuteron is an eigenstate of Sx and not of Sz- We express the 
states \Mx = ±1) in terms of eigenstates of as: 



IM 



X 



±1) = i[|M, = +1) ± V2\Mz = 0) + |M, = -1)]. (35) 



2' 

Then, as before, for A = ±1 the matrix elements for (^)f^ are summed as: 



^dQ' 



J2 \{M',e'\T,NN\e = ex,Mx = +l)\'. (36) 

e',M' 

For (^)-|-^, we simply replace |Mj = by \Mx = —1) in the initial state. 

It is worth mentioning at this point that in Eqs. (j32j) and (|34j) the denominators are not 
the total unpolarized differential cross- sect ion. This is a manifestation of a spin-1 target. 
In the case of a spin-| target they would represent the total unpolarized differential cross- 
section. 
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IV. STRATEGY 



The 7(i elastic scattering process is sensitive only to the isoscalar component of the 
nucleon polarizabilities. For the purpose of this paper, we shall assume that the proton 
polarizabilities are established and focus on the neutron polarizabilities. We believe that 
the debate is not over whether the neutron polarizabilities exist or not, but how accurately 
they can be extracted. Hence, the goal of this project is to build a strategy to facilitate the 
accurate extraction of and jin ■ ■ ■'jAn- The 0{Q^) predictions from Eq. (jH) give a 

starting point in the search. 

An examination of Eq. (j221) reveals that the leading-order term is the Thomson 

term. Since this term is absent for the neutron, its polarizabilities express themselves via 
interference with the proton Thomson term in coherent 'yd scattering. We know that the 
unpolarized differential cross-section (dcs) is proportional to the spin-average of the square 
of the amplitude which makes it apparent that and /3„ occur at O^u"^) and the 7„'s occur 
at in the dcs. This is also true for E. Hence we can say that as we increase the 

photon energy we should first access a„, then /3„ (as it is of order 10% of the value of a„) 
then finally the 7„'s in the dcs or in E. In this paper we focus on S only, as much work has 
already been done on the dcs [l3,[l3l- We might expect to be able to extract an and /5„ at 
intermediate energies (~100 MeV) and focus on the 7„'s at the highest possible energies (in 
our case it is 135 MeV as we want to remain below the pion-production threshold). Since 
we are limited by the energy range, it is only logical to use S to probe the polarizabilities 
it is most sensitive to, namely a„ and However, we shall see that this strategy is not 
guaranteed to work as cancellations occur in S. For the 7„'s we look at other observables, 

and S^, as the situation is more promising there. In S^. and S^, f3n and the spin 
polarizabilities occur at the same order {0{uj^)). Thus, having foreknowledge of a„ and 

and E^; seem ideal for the purpose of extracting the 7„'s. 

We are interested in knowing how accurately the polarizabilities can be extracted and 
at 0{Q^) there are no free parameters. The 0{Q^) predictions for the polarizabilities (see 
Eq. are modified by contributions of higher orders in the chiral expansion 24,|33|- To 
account for the 0{Q^) and higher order effects in a„, and 7i„ . . . 74^, we introduce six 
new parameters as corrections to the 0{Q^) values, which we call Aa„, A/3„ and A7j„(i = 
1,2,3,4). For one particular plot only one of these parameters, for instance, Aa„ is varied 
with the rest being set arbitrarily to zero. This gives us a measure of the sensitivity of that 
particular observable to Aa;„,. 



V. RESULTS AND DISCUSSION 

A. Photon/Beam Polarization Asymmetry 

At 0{Q^) the only term that contributes to E is the Thomson term which has a e' ■ e 
structure. Consequently, the photon/beam polarization asymmetry at leading-order is: 

._(cos'^»-l) (3,) 



(C0S2 ^ + 1) ' 



This result is modified by interference of the Thomson term with the next-to-leading order 
{0{Q^)) term in the -iNN 7A^A^ amplitude. Note that any term in the 7A^A^ 7A^A^ 
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amplitude that has a e ' ■ e structure merely alters the strength of the Thomson term in the 
amplitude and will, on its own, still just lead to the form in Eq. ()37p for the photon beam 
polarization asymmetry, S. 



lA Calculation 




30 60 90 120 150 180 



cm. angle (deg) 



O(Q^) Calculation 




30 60 90 120 150 180 



cm. angle (deg) 

FIG. 5: Plots of S, in Impulse Approximation and to 0{Q^) at 
different energies. 

Fig. El shows the photon polarization asymmetry E, calculated at various energies in 
Impulse Approximation and to 0{Q^) respectively. The difference between these two plots 
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is that the top one does not have any two-body current in the calculation, i.e., T^^n 
arbitrarily set to zero. Comparing the two it is fairly obvious that the two-body currents 
play very little role in this observable up to energies of around 70 MeV. One might wonder 
why the effects of two-body currents are apparently so small. 

As alluded to before, the leading-order result in the chiral expansion) is modi- 

fied by interference with the next-to-leading order {0{Q^)). For example, Eq. (|3ip can be 
expanded as: 



Y: \{M',e'\T%\M,e = l)\'+2^ {M' , e'\T$j,\M, e = 1){M, e = l\T%\M' , e')]+. . . 



(38) 

where the numbers in the superscript denote the chiral order and i = {x,y). The first term 
is the Thomson term which produces the result (jHTj) and the second term is the interference 
between the Thomson term and the next-to-leading order term. Since the Thomson ampli- 
tude is non-zero only for M = M', the interference is sensitive to deuteron spin transitions 
where AM = 0. In contrast, the two-body currents dominate the M = —1 ^ M' = 1 
transition of the deuteron {i.e. AM = 2). But what is the strength of the two-body currents 
for the M = M' spin transition amplitudes for the 'yNN — >■ 'fNN process? To answer this 
question we considered the simplest case involving an S-wave deuteron. Focusing on the 
AM = spin transitions, we found that all the three cases (1 — > 1, — >■ 0, — 1 ^ —1) have 
the same energy behavior and are also roughly of the same strength. 1^1 was chosen as 
the representative because of the relative simplicity in the analytical calculations. To make 
this particularly simple we also took the limit k = k' ^ 0. In this limit diagrams a, c and 
d from Fig. |21 do not contribute to the amplitude at all. Only two of the diagrams in Fig. El 
(b and e) contribute and their contributions are: 

•^{b),{e) = COS 6* (e), 

A^y _ n 
•^(6),{e) - U, 

AV'' — n 

•^(6),(e) - U, 

^S,(e) = ^me)- (39) 

The first (second) letter in the superscript denotes the photon polarization in the initial 
(final) state. ^{b),{e) are two different functions of p,p' and m^r. The result (jHUI) is obviously 
of the form e' ■ e. Of course the approximations used to obtain Eq. (j39|) break down, but 
numerically we have checked that they are reasonable for a; < 70 MeV. Thus, the two-body 
currents do not contribute to S up to photon energies of around 70 MeV. The fact that 
the amplitudes in Eq. ()39|) have a Thomson-like behavior has been computationally verified 
where the same assumptions and simplifications were used. Hence, it is safe to conclude 
that at energies below and up to 70 MeV, an Impulse Approximation calculation does a 
fair job in predicting the photon/beam polarization asymmetry. In practice, however, we 
cannot neglect the two-body currents as they may play a significant role in the differential 
cross-section or when specific spin transitions are considered. Also, our aim is to calculate 
up to the pion-production threshold and the two-body currents are not negligible over the 
full energy range. 

Fig. ini shows the effect of varying a;„ and /?„ on S. In these plots the values of a„ and 
Pn are varied around their central 0{Q^) values, i.e., Aa„ i^Pn) = correspond to the 
value of On (/5n) given in Eq. Q. Varying the polarizabilities in such a manner gives us 
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With varying Aa^ 
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FIG. 6: Plots of S for varying Aa„ and A/?„ at 70 MeV. 

an idea as to how accurately these values can be measured. The first observation is that 
the overall nature of E at 0{Q^) is not very different from what is obtained at (^(Q^), i.e., 
Eq. ()37p . Another obvious observation is that S is not appreciably sensitive to variations 
in either an or The difference in scales of the top and the bottom plot of Fig. E] should 
especially be noted. For < 70 MeV we can use Taylor expansions of the amplitudes (as 
in Eq. (j22j) ) and calculate S in Impulse Approximation in the cm. frame. Then by using 
a Taylor expansion in u for ^ we can analyze the role of a and /3, the isoscalar nucleon 
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polarizabilities [3J|. The result is: 
S I ^ 4: cos 9 



So M (cos2 ^ + 1) 

^2 1 

- TT^TTT ^ TT^l^M^P COS ^(COS^ ^ + 1) + 2(6 COS^ 9 - COS^ 9-7 COS^ ^ - COS ^ + 3) 

2(cos^ 9 + 1)^ 

+ 4/tp(cos2 ^ + 1) (1 - 2 cos 0) - /t5(3 cos^ ^ + 10 cos^ 9- A cos^ + 10 cos ^ - 7) 

+ 4«;3(cos'e + l)(l — COS0 — cos^ 9) + 2k^(cos^ 6* + 1)(1 — cos^^ — cos^ 9) 

+ 4:KpkI{cos'^9 + 1){1-cos9~cos^9) + {kI + kI){1-cos'^9)] + 0{uj^), (40) 

where, Kp{Kn) is the anomalous magnetic moment of the proton (neutron), M is the mass 
of the nucleon, and u is the photon energy. There are three points to be taken away from 
Eq. (gO))- 

1. P occurs at O^u"^) but vanishes at 6' = |, as verified by the bottom plot of Fig. El 
Moreover, according to Eq. (jlj), /3 is approximately one-tenth the value of a and so its 
impact on E is small, even away from |. 

2. a occurs at the next higher order {0{uj^)) and hence is suppressed by an additional 
factor of uj. This may seem peculiar since Eq. ()22|) suggests the occurrence of both a 
and f5 at the same order. The absence of a at 0(a;^) in Eq. (jlUI) is because the terms 
containing a in Eq. ()22|) have an e ' ■ e structure and the a dependence cancels between 
the numerator and the denominator in E to (9((x'^). 

3. Given the second point, one might wonder why the sensitivity to a seems greater than 
p. The next term in the Taylor expansion of ^ contains a but is multiplied by to/M. 
Thus, the ratio of the terms containing a and (3 is proportional to 

_i_ = _._.10x7%~70%. (41) 

Thus, given that the range of variation of Aq;„ is around four times the range of 
variation of A/3„, the sensitivity to appears to be greater. 

The plots in Fig. IHlare a consequence of these facts. The photon/beam polarization asym- 
metry shows very little sensitivity to a„ and /3„ at 70 MeV. Hence, the conclusion is that — 
at least at this energy — S would not be helpful in extracting the values of a„ and with 
the current levels of experimental precision which is ~ 1% at HI7S, especially if we are 
looking to measure a very small variation in a„ from its 0(Q^) predicted value. 

Fig. Elshows the same observable S, but at a higher energy (100 MeV), and its sensitivity 
to and Pn- As expected, the sensitivity grows with energy. However, the spin polariz- 
abilities start playing a role from ~ 90 MeV and it is perhaps not advisable to use E to 
extract and /3„ at these energies. 

The sensitivity to the spin-independent isoscalar polarizabilities is actually much stronger 
in the unpolarized differential cross-section , and so we would advocate using unpolarized 
Compton scattering on the deuteron to determine a„ and 
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FIG. 7: Plots of S for varying Aa^ and A/?^ at 100 MeV. 



B. Double-Polarization Asymmetry 



The spin polarizabilities occur at 0{u^) in the one-body invariant amplitudes. This 
means that they express themselves in the differential cross-section via products like 
^{A*Aj), = 3. . .6, and consequently they first appear at 0{uj^), i.e., at next-to-next- 
to-next-to-next-to-leading order (N^LO) in u. However, when one considers the double- 
polarization asymmetries, the leading-order is of 0{uj^) and not 0{uj'^). (Refer to P] for the 
7P case.) Here, the spin polarizabilities express themselves at 0{uj^) i.e. NNLO, which means 
their effects should be seen at much lower energies than in the unpolarized cross- sect ion. It 
is also true that an and occur at NNLO in the double-polarization asymmetries and one 
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might expect an to have a greater effect than the spin polarizabihties because it is a few 
times larger than the 7n's. However, the 0{uj^) interference terms containing and /3„ are 
suppressed by an additional factor of 1/M (see Eq. ((221)) and the effect of variation in these 
terms is of the same magnitude as from the terms containing the 7„'s, as long as a„ and the 
7„'s are varied by the same amount. Thus, the conclusion is that the double-polarization 
asymmetries are equally sensitive to all of the polarizabihties. This is good news because 
there are other methods by which and Pn can be unambiguously extracted and so at this 
point we can assume that a„ and /3„ have been determined and we are interested in extract- 
ing the spin polarizabihties alone. It is worth noting here that at low energies {uj < 
the deuteron observables are dominated by the proton because of the large contribution 
from the Thomson term. Hence, at these energies, one would expect these observables to 
look similar to the results for the proton. It is indeed the case that at 70 MeV our results 
are similar to those of Ref. for the proton. 



AYj^ (10 fm ) varying Ay,^ (10 fm ) varying 




Lab Angle (deg) Lab Angle (deg) 



FIG. 8: Plots of with varying 7n's at 135 MeV (lab) energy. 



The first two sets of plots (Figs. |H1 and ^ relate to the parallel target polarization asym- 
metry. The plots in Fig. El only show the numerator of Eq. ((321) • It is helpful to look at 
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FIG. 9: Plots of with varying 7„'s at 135 MeV (lab) energy. 
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such a plot because it gives a quantitative estimate of the experimental precision required, 
but, experimentally speaking, it is better to measure as taking the ratio eliminates a 
number of systematic uncertainties. The top left plot results when 71^ is varied, the top 
right corresponds to variation in 72n, the bottom left corresponds to the variation in 73^ and 
the bottom right corresponds to the variation in 74^. For these and all subsequent plots the 
values of the neutron spin polarizabilities have been varied around their OiCf") values, z.e., 
^7m = (z = 1 . . . 4) corresponds to the values in Eq. (^j). 

Except for 7i„, it seems that there is not an appreciable sensitivity to the neutron spin 
polarizabilities in S^. This is more or less obvious because the sensitivity to the spin po- 
larizabilities is at NNLO in uj. The 0((^^^ values of the spin polarizabilities predict that 
7in > 72n > 73n(= " 74n), which suggests that the sensitivities will also be in that order, a 
fact supported by Figs. |H1 and IHl Also, since the range of tSrun (i = 1 . . . 4) chosen here is 
±7m, the value of 71^ was varied the most. However, complications do occur because these 
observables are sensitive to different combinations of the spin polarizabilities. 



Figs. ^1 and ^2 show the perpendicular polarization asymmetry with Fig. ^2 showing 
only the numerator of Eq. (j34j) . Again the neutron spin polarizabilities are varied around 
their O(Q^) values. As before, this observable is not very sensitive to variations in the 
neutron spin polarizabilities, except for 7i„ and 73^. There is a fair amount of sensitivity 
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to 7in, but ignorance of 73„ could hinder extraction of 71^. The projected HI7S error-bar 
for with the current flux is ~ 5%. Thus, with the upgraded flux at HI7S, an experiment 
to extract some combination 71^ and 73^ with good precision from a measurement of S^; 
in elastic scattering seems feasible. This combined with measurements for 70^ and 7,rn 
such as mentioned in Eqs. (fT2|) and |T5|l may make it possible to pin down 7i„, 73^ and a 
combination of 72„ and 74„. 



C. Theoretical Ambiguities 

The theoretical tool used in our calculations is HBxPT in which baryons are introduced 
with the assumption that their rest four momenta dominate their total four momenta, in 
other words, the nucleons are virtually at rest. However, in reality this is not true and 
because of this the pion-production threshold is not in the correct place in the theory. This 
assumption creates a problem for all processes involving baryons near the pion-production 
threshold. For instance, the 7p cross-section varies rapidly near the threshold and if one is to 
believe the HBxPT calculations around this energy, it is crucial to get the pion-production 
threshold in the right place (see, e.g., Ref. |^). The same is true for our calculations too. 
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FIG. 11: Plots of Aa; with varying 7n's at 135 MeV (lab) energy. 
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To remedy this problem Hildebrandt et ah, [3, 



used 



2M' 



(42) 



as the energy 'going into' the amphtude as this includes the nucleon kinetic energy which 
is otherwise neglected in HBxPT. 

To determine the appropriate photon energy to use in Eqs. ()19|) , (jlHI) and (PT|) we examine 
the integrands in the one-body and the two-body currents and extract the energy at which 
they first encounter a pole. To do this, firstly we need to define the kinematics. Working 
in the 7(i CM frame and defining = {u, k) and k'^ = {uj', k') to be the four-momenta of 
the incoming and the outgoing photon, the four-momenta of the outgoing deuteron will be 
, —k) and K'^ = {—B + ^jj-, —k '), where B is the deuteron binding energy. 



Next, 



for the one-body currents, we put the spectator nucleon on its mass shell and calculate 
the four momenta of the struck nucleon using four-momentum conservation. In the 
final step, using these momenta, we extract the energy at which the integrand first 
encounters a pole, which is 

,2 



UJ 



-B + m^ + —f- + 0{m. 
ZMd 



(43) 
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• for the two-body currents, one of the nucleons is put on its mass shell and the four- 
momenta of the nucleon and of the pion(s) are calculated using four-momentum con- 
servation. As in the one-body case, the energy at which the integrand encounters a 
pole is extracted. Interestingly, the energy turns out to be the same as in Eq. 

Hence, for the purpose of our calculations, we redefine 

as the energy 'going into' the amplitude. In all of the results reported in the previous 
two subsections, the value of uj was redefined according to Eq. (jUj). 
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FIG. 12: Plots of at 135 MeV (lab) energy. The dashed 
curve is the HBxPT calculation where u) is simply the photon 
energy in the jd cm. frame. The dot-dashed curve symbolizes 
the other extreme where the nucleons are assumed to be free. The 
solid curve represents the calculation done with Co redefined as in 
Eq. jm). 

Fig. El is a representation of how observables are affected by altering the value of the 
CM energy. The dashed curve is the HBxPT calculation where the nucleon is assumed to 
be heavy and hence practically at rest, i.e. uj = u. The dot-dashed curve symbolizes the 
other extreme where the nucleons are assumed to be free and so oj = u! + The real world 
(for nucleons bound inside a deuteron nucleus) is somewhere in between and the solid curve 
represents such a situation. The uncertainty arising from the freedom to adopt different 
prescriptions for Co is an effect of order 0{Q^) but it obviously has a significant effect on 
our conclusions. The impact of this ambiguity is energy-dependent and largest at 135 MeV. 
For a realistic estimate of this uncertainty, one should consider the discrepancy between the 
solid and the dashed curve. 

Based on NN potentials that describe the low-energy A^A^ data quite accurately, there 
are several deuteron wavefunctions on the market. All of these wavefunctions differ from one 
another in facets like the deuteron D-state probability. Because of these inherent differences, 
the results reported in this paper will be sensitive to the choice of deuteron wavefunctions. 
This too is a higher order effect and should typically enter the result for the jd process 
only at 0{Q^). We have used the wavefunctions obtained from the Nijm93 potential model, 



Free nucleon case 
CO <— cm. energy 




± 



22 



which gives a very good description of the np data to energies of 350 MeV with a P^r 
datum of 1.87 |36|. Fig. CHI gives a flavor of the effect that choosing a different wavefunction 
has on a representative observable. The solid curve uses the next-to-leading order chiral 
wavefunction [s^, whereas the dashed curve uses the Nijmegen wavefunction. It was shown 
in Ref. that these two wavefunctions provide the extrema between which the uncertainty 
due to the choice of wavefunctions lie. It should be noted that the uncertainty is largest for 
the case presented here, i.e., for 135 MeV. 



0.4 



0.2 
" 
-0.2 



"^•"^0 60 120 180 

Lab Angle (deg) 

FIG. 13: This figure shows T^x at 135 MeV (lab) energy for two 
different wavefunctions. 

Comparing Figs. ^1 and ^1 one might be led to the conclusion that though the theory 
predicts some sensitivity to 7i„, the combined theoretical ambiguities may nullify the effect. 
However, it is fair to say that the errors presented above depict the largest uncertainties 
possible and hence give a very conservative estimate of the typical size of the theoretical 
uncertainties. A reasonable range of the error for S^. in the area where there is maximum 
sensitivity to 71^ is ~ 13.8% from the uncertainty because of the choice of wavefunction and 
~ 1.8% due to the ambiguity in the position of the pion-production threshold in HBxPT. 
Combining these two uncertainties in translates into ~ 50% uncertainty in the extraction 
in the value of 7i„. However, much thought is being given toward resolving these issues and 
hence with the passage of time there can only be reduction in the sizes of the theoretical 
ambiguities discussed in this Section. 

VI. CONCLUSION 

From the discussion in the first part of Section |V] we conclude that the photon/beam 
polarization asymmetry is not a good observable to extract a„ and This observable 
is not sensitive to and Pn because effects of a„ cancel in this observable at 0{uj^) and 
sensitivity to /5„ disappears at 9 = ^ where S is largest. However, the second part of 
Section El leads us to believe that the double-polarization observables, especially T,^, can 
be instrumental in extracting information about the neutron polarizabilities. There is real 
hope in this direction, as the HI7S upgrade program will improve the experimental precision, 
while focus on resolving the theoretical ambiguities will reduce the uncertainties discussed 
in Section IV CI 
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However, it is evident that the photon polarization asymmetry and double-polarization 
asymmetries are not sufficient by themselves to extract all the neutron polarizabilities. We 
advocate first trying to extract a„ from measurements of the isoscalar polarizabilities from 
unpolarized cross-section measurements on the deuteron. Then the dispersion sum rule (|iup 
can be used to extract This is not an entirely desirable method because the sum rule 
for the neutron has itself been derived from the deuteron information and contains model 
assumptions. But, it should provide a starting point. 

With the spin polarizabilities the case is worse because the proton spin polarizabilities are 
yet to be established. However, plans to measure spin-dependent polarizabilities using 7p 
scattering are in the pipeline j38|. Hopefully, soon spin polarizability measurements for the 
proton produce promising data. Even with the values for 7ip . . . in hand, measurements 
of 7(i — >• 7(i are not guaranteed to give accurate results, although we have found evidence 
which suggests there will be sensitivity to neutron spin polarizabilities, especially 7i„, in the 
'jd 7(i data. There is also light at the end of the tunnel in the form of 'effective neutron 
targets' like polarized Helium-3, or the use of breakup reactions [3] involving the deuteron 
or Helium-3. Though cross-sections in the breakup reactions are significantly smaller than 
in the elastic case, an enhanced flux at HI7S will surely help the matter. 

To summarize, the long-term plan to extract the neutron polarizabilities should involve 

• extracting a„ and /5„ from unpolarized experiments. 

• measuring 7^^ (i = 1 . . . 4) for the proton. 

• pinning down a linear combination of 7i„ and 73^ from Tj^ measurements on deuteron. 

• further work on 70 and 'Jt^ for both proton and neutron. 
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